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Higgs inflation is a simple and elegant model in which early-universe inflation is 
driven by the Higgs sector of the Standard Model. The Higgs sector can support 
early-universe inflation if it has a large nonminimal coupling to the Ricci spacetime 
curvature scalar. At energies relevant to such an inflationary epoch, the Goldstone 
modes of the Higgs sector remain in the spectrum in renormalizable gauges, and 
hence their effects should be included in the model's dynamics. We analyze the 
multifield dynamics of Higgs inflation and find that the multifield effects damp out 
rapidly after the onset of inflation, because of the gauge symmetry among the scalar 
fields in this model. Predictions from Higgs inflation for observable quantities, such 
as the spectral index of the power spectrum of primordial perturbations, therefore 
revert to their familiar single-field form, in excellent agreement with recent measure- 
ments. The methods we develop here may be applied to any multifield model with 
nonminimal couplings in which the M fields obey an SO{M) symmetry in field space. 
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I. INTRODUCTION 

The recent discovery at CERN of a scalar boson with Higgs-like properties pQ heightens 
the question of whether the Standard Model Higgs sector could have played interesting roles 
in the early universe, at energies well above the electroweak symmetry-breaking scale. In 
particular, the suggestive evidence for the Higgs boson raises the possibility to return to an 
original motivation for cosmological inflation, namely, to realize a phase of early-universe 
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acceleration driven by a scalar field that is part of a well-motivated model from high-energy 
particle physics [2HI]. 

Higgs inflation [5J represents an elegant approach to building a workable inflationary 
model based on realistic ingredients from particle physics. In this model, a large nonminimal 
coupling of the Standard Model electroweak Higgs sector drives a phase of early-universe 
inflation. Such nonminimal couplings are generic: they arise as necessary renormalization 
counterterms for scalar fields in curved spacetime [6H9]. Moreover, renormalization-group 
analyses indicate that for models with matter akin to the Standard Model, the nonminimal 
coupling, £, should grow without bound with increasing energy scale [9]. Previous analyses of 
Higgs inflation have found that £ typically grows by at least an order of magnitude between 
the electroweak symmetry-breaking scale and the inflationary scale [TUHT2] . 

The Standard Model Higgs sector includes four scalar degrees of freedom: the (real) Higgs 
scalar and three Goldstone modes. In renormalizable gauges, all four scalar fields remain 
in the spectrum at high energies [T2TH6] . Thus the dynamics of Higgs inflation should be 
studied as a multifield model with nonminimal couplings. An important feature of multifield 
models, which is absent in single-field models, is that the fields' trajectories can turn within 
field space as the system evolves. Such turns are a necessary (but not sufficient) condition 
for multifield models to depart from the empirical predictions of simple single-field models 



In this paper we analyze the background dynamics of Higgs inflation, in which all four 
scalar fields of the Standard Model electroweak Higgs sector have nonminimal couplings. 
We find that multifield dynamics damp out quickly after the onset of inflation, before per- 
turbations on cosmologically relevant length scales first cross the Hubble radius. As regards 
observable quantities like the power spectrum of primordial perturbations, the model there- 
fore behaves effectively as a single-field model. The multifield dynamics remain subdominant 
in Higgs inflation because of the particular symmetries of the Higgs sector. Closely related 
models, which lack those symmetries, can produce conspicuous departures from the single- 
field case 123]. 



We are principally interested here in the behavior of classical background fields and 
long-wavelength perturbations, which behave essentially classically. Therefore we bracket, 
for this analysis, the question of the unitarity of Higgs inflation. Conflicting conclusions 
have been advanced regarding whether the appropriate renormalization cut-off scale for this 
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model should be M p \, M p i/^, or M p i/£, where M p \ = (SixG)^/ 2 is the reduced Planck mass 
[L2l HU [151 El] . Even if Higgs inflation might conclusively be shown to violate unitarity, the 
techniques developed here for the analysis of multifield dynamics will be relevant for related 
models that incorporate multiple scalar fields with nonminimal couplings and symmetries 
(such as gauge symmetries) that enforce specific relations among the couplings of the model. 
In particular, we expect that multifield effects in models with M scalar fields, in which the 
scalar fields obey an SO(N) symmetry, should damp out rapidly. 

In Section II, we briefly introduce the multifield formalism and establish notation. We 
apply the formalism to Higgs inflation in Section III, and in Section IV we analyze the 
behavior of the turn-rate, which quantifies the rate at which the background trajectory of 
the system deviates from a single-field case. We study how quickly the turn-rate damps 
to zero, both analytically and numerically, confirming that for Higgs inflation the turn-rate 
becomes negligible within a few efolds after the start of inflation. In Section V we turn 
to implications for observable features of the primordial power spectrum, confirming that 
multifield Higgs inflation reproduces the empirical predictions of previous single-field studies. 
Concluding remarks follow in Section VI. 



II. MULTIFIELD DYNAMICS 

Following the approach established in [23], we consider models with M scalar fields in 
(3+1) spacetime dimensions. We use Greek letters to label spacetime indices, /i, v — 0, 1, 2, 3; 
lower-case Latin letters to label spatial indices, i,j = 1,2,3; and upper-case Latin letters 
to label field-space indices, J, J = 1, 2, M . We also work in terms of the reduced Planck 
mass, M p i = (87rG) -1 / 2 . In the Jordan frame, the action takes the form 



^Jordan — / gT^Ta/— f 



f(tf)R - Uug^d^d^ 3 - V& 1 ) 



(1) 



Here is the nonminimal coupling function, and we use tildes for quantities in the 

Jordan frame. We perform a conformal transformation to the Einstein frame by rescaling 
the spacetime metric tensor, 

2 

9^( x ) = Jpf^\ x )) 9^(x), (2) 
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so that the action in the Einstein frame becomes [25 

' M 2 



Einstein 



d Xy 



-g 



1 



(3) 



The potential in the Einstein frame, V, is related to the Jordan-frame potential, V, as 



4/ 2 



(4) 



and the coefficients of the noncanonical kinetic terms are (25J |26] 



Ml, 



fir) 



(5) 



where fj = df /dcj) 1 . The nonminimal couplings induce a field- space manifold in the Einstein 
frame that is not conformal to flat; Qjj serves as a metric on the curved manifold [25J. 
Therefore we adopt the covariant approach of [23], which respects the curvature of the 
field-space manifold. 

Varying Eq. (p| with respect to / yields the equation of motion, 



□0 7 + gvr I JK d ll <f> J d v <f - G IK V K = 0, 



(6) 



where D0 7 = g^^^ y and T I JK ((f) L ) is the Christoffel symbold for the field-space manifold, 
calculated in terms of Qjj. We expand each scalar field to first order around its classical 
background value, 

0V) = /W + ¥V), (7) 

and also expand the scalar degrees of freedom of the spacetime metric to first order around 
a spatially flat Friedmann-Robertson- Walker (FRW) metric [27H29] 



ds 2 = -(1 + 2A)dt 2 + 2a{d i B)dx i dt + a 2 [(1 - 2VO<& - + 2d i d j E] dx*dx j , 



(8) 



where a(t) is the scale factor. We further introduce a covariant derivative with respect to 
the field-space metric and a directional derivative along the background fields' trajectory, 
such that for any vector A 1 in the field-space manifold we have 



VjA 1 = d J A I + T I JK A K , 

V t A' = tfVjA 1 = A 1 + r T JK A J ip K , 

where overdots denote derivatives with respect to cosmic time, t. 



(9) 



5 



To background order, Eq. (|6| becomes 



v^ 1 + m^ 1 + g IK v K = o, 



(10) 



projects onto the subspace orthogonal to a 1 . Eq. (10) then simplifies to 



where all quantities involving Gu(<j) K ), V^ 1 ), and their derivatives are evaluated at back- 
ground order in the fields: Qu — > Qu(f K ) and V — > V(ip T ). Following [T7] we distinguish 
between adiabatic and entropic directions in field space by introducing a unit vector 

(id 

where 

a = = yjGuy 1 ^. (12) 

The operator 

g" = gu _ (13) 



a + 3H& + V a = 



(14) 



where 

V a EE (7%. 

The background dynamics likewise take the simple form 



H 2 
H 



3M p 2 



-a 2 + V 



2M p f 2 ' 



(15) 



(16) 



where H = a/a is the Hubble parameter. 

We may also separate the perturbations into adiabatic and entropic directions. Working 
to first order in perturbations, we introduce the gauge-invariant Mukhanov-Sasaki variables 
[27H30] 

Q 1 ee Stf + £j> (17) 

and the projections 



5s 1 ee s^Q- 1 . 



(18) 
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The gauge-invariant curvature perturbation may be defined as 1Z C 
where the perturbed energy-momentum flux is given by T° i = di5q 
that 1Z C is proportional to Q a [23J 



a 



4 - [ H /(p + p)]Sq, 
We then find 

(19) 



Expanding Eq. (|6j) to first order and using the projections of Eq. (18), the perturbations 
Q a and 5s 1 obey 



Q a + 3H Q a + 



+ M C 



d 



u 2 - 



1 d ( a 3 & 2 



M p 2a 3 dt V H 



2 Jt M sJ ) 



2|^ 



H 
H 



(coj5s J ) 



(20) 



and 



V 2 5s T + [3115^ + 2& I coj] Vtbs 1 + 



where the mass-squared matrix is 



%5 I J + M I j-2a I (M aJ + -ujj 



Qa T tQ& • Q (J 

H a 



5s" 



(21) 



m'j = G IK {VjV K v) - n'jMjiffvF, 

M aJ = aiM 1 ], M aa = a I a J M I J . 



The turn-rate [221 [23] is given by 



L) 



Vtcr 1 



-Vks IK , 
a 



(22) 



(23) 



and uj = |a; J |. Eqs. (20) and (21) decouple if the turn-rate vanishes, u 1 = 0. In that case, 



Q a evolves just as in the single-field case 



. Given Eq. (19), that means that 



the power spectrum of primordial perturbations, Viz, would also evolve as in single-field 
models. Thus a necessary (but not sufficient) condition for multifield models of this form 
to deviate from the empirical predictions of simple single-field models is for the turn-rate to 
be nonnegligible for some duration of the fields' evolution, u 1 ^ 0. 



7 



III. APPLICATION TO HIGGS INFLATION 

The matter contribution to Higgs inflation [S] consists of the Standard Model electroweak 
Higgs sector, which may be written as a doublet of complex scalar fields, 



h= (24) 




The complex fields h + and h° may be further decomposed into (real) scalar degrees of 
freedom, 

f (25) 
h = — (4> + iX 3 ), 

where <fi is the Higgs scalar and \ a (with a = 1, 2, 3) are the Goldstone modes. In the Jordan 
frame, the potential V^ 1 ) depends only on the combination 

tfh = \[4> 2 + X 2 } , (26) 

where x — (xSx^X 3 ) is a 3-vector of the Goldstone fields. In particular, the symmetry- 
breaking potential may be written 

V(<P I ) = ^(<P 2 + X 2 -v 2 ) 2 , (27) 

in terms of the vacuum expectation value, v. For the Standard Model, v = 246 GeV <C M p \. 
For Higgs inflation, the nonminimal coupling function is given by 

m = ^ + &h = I [Ml + £ (0 2 + x 2 )] , (28) 

where Mq = M 2 X — £v 2 and £ > is the dimensionless nonminimal coupling constant. In 
Higgs inflation, we take £ ~ (9(10 4 ) [5], and therefore we may safely set M 2 = M 2 V In 
the Einstein frame, the potential gets stretched by the nonminimcal coupling function /((j) 1 ) 
according to Eq. Q. Given Eqs. (27) and (28), this yields 



4 ( ±2 i _,2 „,2\ 2 



AM* (0 2 + x 2 ~v 



V^ 1 ) = — plvy — (29) 

4[M 2 +£(^ + x 2V 2 



The model is thus symmetric under rotations among and \ a that preserve the magnitude 



\/(j) 2 + X 2 - When written in the "Cartesian" field-space basis of Eq. (25), in other words, 
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the SU(2) electroweak gauge symmetry manifests as an 5*0(4) spherical symmetry in field 
space. 

For any model with M real- valued scalar fields that respects an SO(Af) symmetry, the 
background dynamics depend on just three initial conditions: the initial magnitude and ini- 
tial velocity along the radial direction in field space, and the initial velocity perpendicular to 
the radial direction. Without loss of generality, therefore, we may analyze the background 
dynamics of Higgs inflation in terms of just two real-valued scalar fields, an d Xi an d we 
may set x(0) = 0, specifying only initial values for 0(0), 0(0), and x(0). This reduction in 
the effective number of degrees of freedom stems entirely from the gauge symmetry of the 
Standard Model electroweak sector. The remaining dependence on x, meanwhile, distin- 
guishes the background dynamics from a genuinely single-field model, in which one neglects 
the Goldstone fields altogether. For the remainder of this paper, we exploit the gauge sym- 
metry to consider only a single Goldstone mode, x ~ >* Xi reducing the problem to that 
of a two-field model. Then /(0 7 ) and V(0 7 ) depend on the background fields only in the 
combination 

r = vV + X 2 - (30) 

Previous analyses [SJ [2SI E2H3I] which considered single-field versions of this model (ne- 
glecting the Goldstone modes) found successful inflation for field values £0 2 ^> M^. We 
confirm this below for the multifield case including the Goldstone modes. The reason is easy 

f 2 
'pi' 



to see from Eq. (29). In the limit £(0 2 + x 2 ) — ^ ^pi? the potential in the Einstein 



frame becomes very flat, approaching 



AM 



4 r / 

pi 



1 + 01 P 1 



ir 2 



(31) 



4£ 2 
See Fig. g) 

Given £ ~ 10 4 , the initial energy density for this model lies well below the Planck scale, 
p ~ V ~ AMpj/£ 2 ~ 10~ 9 Mp[. In fact, as we will see, successful slow-roll inflation (pro- 
ducing at least 70 efolds of inflation) occurs for initial values of the fields below the Planck 
scale, unlike in models of chaotic inflation with polynomial potentials that lack nonminimal 
couplings. Moreover, as emphasized in [5], the flattening of the potential in the Einstein 
frame at large field values makes Higgs inflation easily compatible with the latest observa- 
tions of the spectral index, n s . Ordinary chaotic inflation with a A0 4 potential and minimal 
coupling, on the other hand, yields a spectral index outside the 95% confidence interval for 
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FIG. 1: The potential for Higgs inflation in the Einstein frame, V((p,x)- Note the flattening of 
the potential for large field values, which is quite distinct from the behavior of the Jordan- frame 

potential, V(</),x) in Eq. (I27I). 



the best-fit value of n s |31j . Below we confirm this behavior for Higgs inflation even when 
the Goldstone degrees of freedom are included. 

The field-space metric Qjj is determined by the nonminimal coupling function, 
and its derivatives. Explicit expressions for the components of Qu for a two-field model 
with arbitrary couplings, £^ and £ x , are given in the Appendix of [23] • In the case of Higgs 
inflation, the SU(2) gauge symmetry enforces £^ = £ x = £. Given this symmetry, the 
convenient combination, C(0 7 ), introduced in the Appendix of (23] becomes 



Ctf) =2f + 6e (0 2 + X 2 ) = Ml + £(1 + 6£)r 2 . 



(32) 



For £0 = £ x = £, the Ricci curvature scalar for the field-space manifold, as calculated in 
takes the form 



(33) 



During inflation, when £r 2 3> M^, this reduces to 



(34) 



indicating that the field-space manifold has a spherical symmetry with radius of curvature 
r c ~ y/£ r. As shown in [23], the curvature of the field-space manifold remains negligible in 
such models until the fields satisfy £r 2 <C M^, near the end of inflation. 
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From Eq. (10), and using the expressions for Q IJ and r JK in the Appendix of 
equation of motion for the background field 4>(t) takes the form 

£(i + 60 



+ 3H(f) 



C 



X 



/ + AM pl — 



0. 



the 



(35) 



The equation for \ follows upon replacing 
velocity vector becomes 

■2 



a 



Ml 



» X- Using Eq. ( |12| ), the square of the fields 
3/ 2 



+ 



./ 



and the gradient of the potential in the direction a 1 becomes 



AM* (0 2 + X 2 ) / 



my 



a 



(36) 



(37) 



We may verify that multifield Higgs inflation exhibits slow-roll behavior for typical choices 
of couplings and initial conditions. First consider the single-field case, in which we set 



X = X = 0- Near the start of inflation (with £0 2 3> M 2 j), the terms in Eq. (35) that stem 



from the field's noncanonical kinetic term take the form 

fll + efl ^ /, 2 

The usual slow-roll requirement for single-field models, \<j>\ <C \H<f>\, ensures that the terms 



(38) 



in Eq. (38) remain much less than the 3H(f> term in Eq. (35). Neglecting 0, the single-field, 



slow-roll limit of Eq. ( 35 ) becomes 



3H<f) 



or, upon using H 2 



V/{3M 2 



pill 



6£ 3 </> ' 



AM* 



3^3 e 



(39) 



(40) 



Setting £ = 10 and fixing the initial field velocity by Eq. (40) requires 0(0) > 0.1M p i to 
yield N > 70 efolds of inflation in the single-field case. 

A much broader range of initial conditions yields N > 70 efolds of inflation in the two-field 



case. From Eq. (16) we see that inflation (with a > 0) requires a 2 <C V. Given the SO(J\f) 
symmetry of the model, we may set x(0) = without loss of generality, and parameterize 
the fields' initial velocities as 

h X ' 

(41) 



0(0) 



X(Q) = — ^ y 
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in terms of dimensionless constants x and y. (The single-field case corresponds to x 
—l,y = 0.) Near the start of inflation, when £r 2 = £<f) 2 3> M 2 j, Eq. (36) becomes 

r4x / M 2 ^ - ! 



(j 2 | x (o)=o 



AM, 



pi 



pi 



[(l + 60x 2 + y 2 ]. 



(42) 



4e ; wmj 27^ 

The first term in parentheses is just the value of the potential, V, near the start of inflation 



as given in Eq. (31 ). The second term in parentheses is small near the beginning of inflation, 



given £r 2 3> M 2 V Hence the initial values for <fi and Xi parameterized by the coefficients x 
and y, may be substantially larger than in the single-field case while still keeping a 2 <C V. 
Fig. p\ shows H(t), (f>(t), and x{t) f° r a scenario in which 0(0) and x(0) greatly exceed 



the single-field relation of Eq. (40): \x\ = 10 and \y\ = 10 . As is evident in the figure, the 



large initial velocities cause the fields to oscillate rapidly. The extra kinetic energy makes 
the initial value of H(t) larger than in the corresponding single-field case. The increase in 
H, in turn, causes the fields' velocities to damp out even more quickly, due to the 3H<f> and 
3Hx Hubble-drag terms in each field's equation of motion. Thus the system rapidly settles 
into a slow-roll regime that continues for 70 efolds of inflation. As shown in Fig. 3, we may 
achieve iV > 70 efolds of inflation with even smaller initial field values by making the initial 
field velocities correspondingly larger. 

IV. TURN RATE 



The components of the turn-rate, u 1 in Eq. (23), take the form 



AM* r 2 



& 2f 



C Af 2 a 2 



+ XX 



(43) 



The other component, u x , follows upon replacing < — > x- The length of the turn-rate 
vector is given by 

' ' (44) 



U) = \U) 



a 



s KM V K V 



KV,M 

I iv 



where the final expression follows upon using the definition of u in Eq. (23 ) and the identity 



-KM _ zK ~MA 
i — 6 aS 



which follows from Eq. (13). We find 



•2, ,2 t,KM\r t/ 

a u) = s V K v m 



\ 2 M™ 
J2jfc r ' 



[c - £V] - (V/ CT ) 2 . 



(45) 



The evolution of the turn rate for typical initial conditions is shown in Fig. |4| 
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FIG. 2: The evolution of H(t) (black dashed line) and the fields 0(i) (red solid line) and x(t) 
(blue dotted line). The fields are measured in units of M p \ and we use the dimensionless time 
variable r = vA M p \t. We have plotted 10 3 H so that its scale is commensurate with the 
magnitude of the fields. The Hubble parameter begins large, H(0) = 8.1 x 10~ 4 , but quickly falls 
by a factor of 30 as it settles to its slow- roll value of H = 2.8 x 10~ 5 . Inflation proceeds for 
At = 2.5 x 10 6 to yield N = 70.7 efolds of inflation. The solutions shown here correspond to 
f = 10 4 , 0(0) = 0.1, xM = 0, 0(0) = -2 x 10~ 6 , and x(0) = 2 x 10" 2 . For the same value of 
0(0), Eq. (40) corresponds to 0(0) = — 2 x 10~ 8 for the single-field case. 



In order to analyze the evolution of the background fields, it is easier to move from 
Cartesian to polar coordinates, in which the angular velocity and turn rate have more 
intuitive behavior. In addition to the radius, r 2 = 2 + x 2 , we also define the angle 

X ' 



- = arctai) ( — ) ■ (46) 
Single-field trajectories correspond to constant j(t). In the polar coordinate system, the 



background dynamics of Eq. (|16|) may be written 

1 



H 1 



H = - 



12/ 
1 



r 



\ M 2 

•2 , 2-2 i U S> 2-2, pi 



4/ 



■;i ..2-2 i 3£ 2 2-2 



_|_ ^..jL^a _|_ — - — y- 



(47) 



The equations of motion become 



f + 3Hr — r7 2 + 



c 



r (f 2 + r 2 7 2 ) - V 2 r + XM ] 



(48) 



and 



M 2 
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FIG. 3: Contour plots showing the number of efolds of inflation as one varies the fields' initial 
conditions, keeping £ = 10 4 fixed. In each panel, the vertical axis is x(0) and the horizontal axis 
is 0(0). The panels correspond to 0(0) = 10" 1 M pi (top left), 10~ 2 M pi (top right), 5 x 10~ 3 M pl 
(bottom left) and 10 -4 M p \ (bottom right), and we again use dimensionless time r = yXMpit. In 
each panel, the line for N = 70 efolds is shown in bold. Note how large these initial velocities are 



compared to the single- field expectation of Eq. ( 40 ) . 



In this new basis the turn rate may be written compactly as 



2 



4 • 

r 7 



2fC I r 2 7 2 (M p 2 ! + £r 2 ) + r 2 C 



(50) 



This expression vanishes in both the limits \ j\ — > and |7| — > 00: if the angular velocity is 
either too large or too small, the fields' evolution reverts to effectively single-field behavior 
(either purely radial motion or purely angular motion). Of the two limits, however, only 
pure-radial motion is stable. It is ultimately the evolution of j(t) that will determine the 
fate of the turn rate. 



It is obvious from Eq. (49) that the line 7 = is the fixed point of the angular motion. 



The character of the fixed point is defined by the sign of the 7 term, which is less trivial. 
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x x 10 5 xx 10 4 



FIG. 4: Evolution of the turn rate. The left picture shows the evolution with initial conditions as 
in Fig. [i| The right figure has initial conditions (f)(0) = 0.1, x(0) = 0(0) = 0, and x(0) = 2 x 10~ 5 
in units of M p \ and r = y/XM p it. In both cases we set £ = 10 4 . Recall from Fig. [ijthat inflation 
lasts until r en d ~ O(10 6 ) for these initial conditions; hence we find that uj damps out within a few 

efolds after the start of inflation. 



It can be negative close to r = due to the high curvature of the field manifold and the 
small value of the Hubble parameter, but in the slow-roll regime of the radial field, with 
ir 2 > Mp 1; the sign of 7 is safely positive. That means that we can treat the angular motion 
as damped throughout inflation. 

For large nonminimal coupling and / or slow rolling of the radial field the last term in Eq. 



(49) may be neglected, which yields 

7 + 3#7 = 0. (51) 



The only complicated object in Eq. (51) is the Hubble parameter, which may be simplified 



in the limit of a slow rolling radial field and large nonminimal coupling upon making use of 



Eq. (47): 




. AM n 2 , 
H ^—\H 2 + ^f- ( 52 ) 



Then Eq. (51) becomes 




\M 2 

7 2 + ^0-0. (53) 



Although Eq. (53) can be solved exactly (see the Appendix), it is instructive to examine 



the two limits of large and small 7, which provide most of the relevant information. 
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For small angular velocity, 7 <C we recover the linear limit 



with the solution 



7 + 



7 = 70 exp 



— — 7 ~ 
12 ' 



3A 



M pl t 



oc e 



-3iV 



(54) 



(55) 



where N = Ht. It is very easy to measure time in efolds in this limit, since the Hubble 



parameter is nearly constant. Eq. (55) illustrates that any small, initial angular velocity 



will be suppressed within a couple of efolds, or equivalently within a time of the order of 

e/(VXM pl ). 



In the opposite limit, 7 ^> y which we call the nonlinear regime, Eq. (53) 
becomes 



7 + 3 



7 







with the solution 



7 



1 3t 

,7o V^, 



(56) 



(57) 



Given Eqs. (55) and (57), we may follow the evolution of any initial angular velocity. If 7 



begins large enough it will start in the nonlinear regime, where it will stay until it becomes 
of order JAMV^. We parameterize the cross-over regime as 



7 



AM 



pi 



(58) 



where z is some constant of order one. The cross-over time may then be estimated by 



inverting Eq. (57) to find 



ni 



2£ 1 

A Mpi z 



1 

7o 



(59) 



There exists an upper limit on the time it takes for the angular velocity to decay, namely 

\^M pl t nl , m£lx = ~. (60) 



We have verified all of these analytic predictions using numerical calculations of the 
exact equations for the coupled two-field system in an expanding universe. In Fig. [5] we 
plot the number of efolds from the beginning of inflation at which the turn rate reaches its 
maximum value, as we vary the fields' initial velocities. Note that for any combination of 
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FIG. 5: Contour plots showing the number of efolds at which the maximum of the turn rate 
occurs, as one varies the fields' initial conditions. In each panel, the vertical axis is x(0) and the 
horizontal axis is 4>(0). The panels correspond to 0(0) = KT 1 M pl (top left), 1CT 2 M pl (top right), 
5 x 1CT 3 Mpi (bottom left) and 10~ 4 M p i (bottom right). We set £ = 10 4 and use the 
dimensionless time- variable r = yf\ M p \t. The thick back curve is the contour line of initial 

conditions that yield 70 efolds of inflation. 



initial conditions that yields at least iVtot = 70 efolds of inflation, u reaches its maximum 
value between iV(co> max ) = 3.5 and 5 efolds from the start of inflation (for the range of 
initial conditions considered there). In Fig. [6] we plot w as a function of time from the 
start of inflation, as we vary the initial angular velocity, 7(0)- The curves in red correspond 
to initial conditions in the linear regime, while the curves in blue start in the nonlinear 
regime. Note that the curves starting in the nonlinear regime have the same amplitude. 
The existence of a maximum time, t n i,max, is evident from the bunching of the blue curves. 
We find y/\M pl t n i max — Tni,max ~ few x £ ~ 10 4 , as expected from Ecj. (|60|) . In these units 
and for the initial conditions used in Fig. [6j inflation lasts until r cn d ~ (9(10 6 ), so r n i jmax 
occurs very early after the onset of inflation. 
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FIG. 6: The turn rate as a function of time for different values of the initial angular velocity. 
The parameters used are £ = 10 4 , 0(0) = 0.1M pl , 0(0) = x(0) = 0, and ^|| < 7(0) < in 

terms of dimensionless time, r = \f\Mp\t. In these units and for 0(0) = 0.1M p i, inflation lasts 

Until Tend ~ 0(1O 6 ). 



Eq. (55) shows that the linear region lasts at most a few efolds, so the duration of 
the nonlinear region is what will ultimately determine whether or not multifield effects will 
persist until observationally relevant length scales first cross the Hubble radius. In the 



nonlinear regime, Eq. (52) yields H ~ 7/a/6£ with 7 given by Eq. (57). The number of 



efolds for which the nonlinear regime persists is given by 

N nl = r Hdt ~ -L r 7* = J H \ / ~T T7~~7 1 ^ (6-1 > 







6£ Jo 




We examine Eq. (61) numerically by fixing £ = 10 4 and 0(0) = 0.1M p i and choosing pairs 
of initial velocities, 0(0) and x(0), that yield 70 efolds of inflation (see Fig. ^ left); and 
also by setting 0(0) to various constant values and varying x(0) (Fig. [7J right). The results 
fall neatly along a least-squares logarithmic fit, as expected from Eq. ( |6~Tj ). The function 
N D \ grows slowly. In order for multifield effects to remain important more than a few efolds 
after the start of inflation, the initial angular velocity would need to be enormous: at least 
ten orders of magnitude larger than typical values of the initial field velocity for single- 



field inflation, as given in Eq. (40). We do not know of any realistic mechanism that could 
generate initial field velocities so large. Moreover, for many combinations of initial conditions 
shown in the righthand side of Fig. [7j N tot > 70 efolds (several sets of initial conditions 
yield iVtot ~ 90 efolds). For those scenarios, the turn rate reaches its maximum value deep 
within the early phase of inflation, long before observationally relevant perturbations first 
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FIG. 7: Number of efolds after the start of inflation until the maximum value of the turn rate is 
reached, as a function of x(0). On the left we plot N(uj max ) for initial conditions giving 70 efolds 
of inflation; on the right we plot the same quantity for various values of 0(0). The logarithmic fit 



is an excellent match to our analytic result, Eq. (61). 



cross the Hubble radius. 

As a final test of our analysis we set £ = 10 2 instead of £ = 10 4 . The smaller value of 
the nonminimal coupling does not lead to a viable model of Higgs inflation — the WMAP 
normalization of the power spectrum requires a larger value of £ [S] — but we may nonetheless 
study the dynamics of such a model. We collect the important information about the 
dynamics of this model in Fig. [SJ As expected, the model can provide 70 or more efolds of 
inflation for a wide range of parameters, and the corresponding turn rate peaks well before 
observationally relevant length scales first crossed the Hubble radius, even when we increase 
x(0) to a few hundred in units of r = y/XM p \t. The excellent logarithmic fit of the time 
at which the turn rate is maximum versus x(0) is again evident. Finally the curves of the 
turn rate versus time show the same qualitative and quantitative characteristics as Fig. [6] 
for £ = 10 4 . Specifically, if one rescales time and the turn rate by £, the two sets of curves 
would be hardly distinguishable. 



V. IMPLICATIONS FOR THE PRIMORDIAL SPECTRUM 

We have found that in models with an SO (AT) symmetry among the scalar fields, the turn 
rate quickly damps to negligible magnitude within a few efolds after the start of inflation. 
In this section we confirm that such behavior yields empirical predictions for observable 
quantities like the primordial power spectrum of perturbations that reproduce expectations 
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FIG. 8: Dynamics of our two-field model with £ = 10 2 , </>(0) = lM p i, and x(0) = 0. Clockwise 
from top left: (1) Contour plot showing the number of efolds as one varies the fields' initial 

conditions. The thick curve corresponds to 70 efolds. (2) Contour plot showing the number of 

efolds at which the maximum of the turn rate occurs, as one varies the fields' initial conditions 
The thick curve corresponds to 70 efolds of inflation. (3) Number of efolds after the start of 
inflation until the maximum value of the turn rate is reached for initial conditions giving 70 
e-folds of inflation, along with a logarithmic fit. (4) The turn rate as a function of time for 

different values of the initial angular velocity, with <^>(0) 

'X Mpit. 



and %i < 7(0) < M, in units of 



from corresponding single-field models. 

For models that behave effectively as two-field models, which includes the class of SO(J\f)- 
symmetric models we investigate here, we may distinguish two scalar perturbations: the 



perturbations in the adiabatic direction, Q a defined in Eq. (18), and a scalar entropic 
perturbation [23] . 

Q s = —5s 1 . (62) 



We noted in Eq. (19) that Q a is proportional to the gauge-invariant curvature perturbation, 
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1Z C . We adopt a similar normalization for the entropy perturbation, 

S = -Q s . 
a 

In the long- wavelength limit, the adiabatic and entropic perturbations obey [231 135] 

/ p 

TZ C = aHS + 
S = (3HS + O 
so that we may define the transfer functions 



(63) 



a 2 H 2 

k 2 
a 2 H 2 



(64) 



Tns(t*,t) = J dt' a(t)H(f)T S s{U,f) 



Tss(t*,t) = exp 



dt' P(t')H(t') 



t* 



(65) 



We take to be the time when a fiducial scale of interest first crossed the Hubble radius 
during inflation, defined by a 2 (t >t: )i/ 2 (t*) = k 2 . In [23], we calculated 

2w(t) 



a(t) 



H(t) 



and 



(3(t) = -2e - i] ss + rjaa - 



3^' 



(66) 



(67) 



where e = —H/H 2 and the other slow- roll parameters are defined as 



Vss = M. 



Pi y 

2 ''' 



<jji<jj J M 1 
pl u 2 V 



(68) 



The dimensionless power spectrum is given by 

k 3 



and hence, from Eqs. (64) and (65) 



Vn(k)=Pn(k*) [l + Tl s (U,t 



(69) 



(70) 



where k corresponds to a length scale that crossed the Hubble radius at some time t > t* 
The spectral index is then given by 



n s {t) = n,(U) - [a(t) + f3(t)T ns (t,Q] sin(2A), 



(71) 
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where 

cosA = — =H=. (72) 
In the limit (co/H) <C 77 CT(J , the spectral index evaluated at N* assumes the single-field form 

[2H1 E21 EU, 

n s (Q = l-6e(Q+2rj afT (Q. (73) 

Crucial to note is that the turn rate, u, serves as a window function within T-jis(t,t*): 
once the coefficient a = 2uj/H becomes negligible, there will effectively be no transfer of 
power from the entropic to the adiabatic perturbations, much as we had found by examining 



the source terms on the righthand sides of Eqs. (20) and (21). The question then becomes 
whether u(t), and hence T-jis(t*,t), can depart appreciably from zero at times when pertur- 
bations on length scales of observational interest first cross the Hubble radius. 

The longest length scales of interest are often taken to be those that first crossed the 
Hubble radius iV* = 55 ± 5 efolds before the end of inflation [27T429] . Closer analysis 
suggests that length scales that first crossed the Hubble radius N* = 62 — 63 efolds before 
the end of inflation correspond to the size of the present horizon [36J. Meanwhile, we follow 
[2"S] in assuming that successful inflation requires N tot > 70 efolds to solve the horizon and 
flatness problems. The question then becomes whether u(t), and hence Tns(t*, t), can differ 
appreciably from zero for iV* < 63. Given the analysis in Section IV, the best chance for this 
to occur is for initial conditions that produce the minimum amount of inflation, iV to t — 70. 

In Table I, we present numerical results for key measures of multifield dynamics. In each 
case we set £ = 10 4 , 0(0) = 0.1M p i, and x(0) = 0. We vary x(0) as shown and adjust 0(0) 
in each case so as to produce exactly N tot = 70 efolds of inflation. Because Tus remains so 
small in each of these cases, there is no discernible running of the spectral index within the 
window = 63 to iV* = 40 efolds before the end of inflation. If we consider a fiducial scale 
k* that first crosses the Hubble radius at = 63 efolds before the end of inflation, then we 
find n s = 0.97 across the whole range of initial conditions, in excellent agreement with the 
measured value of n s = 0.968 ± 0.012 [31]. If instead we set k* as the scale that first crossed 
the Hubble radius N* = 60 efolds before the end of inflation, we find n s = 0.967 across the 
entire range of initial conditions. 
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x(o) 


w(iV* = 63) 


T-jis (max) 


n s (N* = 63) 


n 8 (N* = 60) 


io- 2 


1.16 x 1(T 10 


2.68 x 10~ 6 


0.969 


0.967 


10- 1 


1.20 x 10" 9 


2.76 x 10~ 5 


0.969 


0.967 


1 


9.41 x 10~ 9 


2.18 x IO" 4 


0.969 


0.967 


10 1 


1.18 x 10~ 7 


2.72 x 10~ 3 


0.969 


0.967 


10 2 


1.12 x 10~ 6 


2.59 x 10~ 2 


0.973 


0.967 



TABLE I: Numerical results for measures of multifield dynamics for Higgs inflation with £ = 10 4 . 

We use dimensionless time r = \fXM v \t. 

VI. CONCLUSIONS 

In this paper we have analyzed Higgs inflation as a multifield model with nonminimal 
couplings. Because the Goldstone modes of the Standard Model electroweak Higgs sector 
remain in the spectrum at high energies in renormalizable gauges, we have incorporated 
their effects in the dynamics of the model. Because of the high symmetry of the Higgs 
sector — guaranteed by the SU(2) electroweak gauge symmetry, which manifests as an 
SO (A) symmetry among the scalar fields of the Higgs sector — the nonmiminal couplings 
for the various scalar fields take precisely the same value = £ x = £) , as do the tree- level 
couplings in the Jordan-frame potential (A^ = A x = A, and so on). The effective potential in 
the Einstein frame therefore contains none of the irregular features, such as bumps or ridges, 
that were highlighted in [23] for the case of multiple fields with arbitrary couplings. With 
no features such as ridges off of which the fields may fall during their evolution, Hubble 
drag will always cause any initial angular motion within field space to damp out rapidly. 
Increasing the initial angular velocity to arbitrarily large values — well into what we call the 
nonlinear regime — only increases the value of H at early times, which makes the Hubble 
drag even more effective and hence hastens the damping out of the multifield effects. 

The rapidity with which the turn-rate damps to zero combined with the requirement 
of N tot > 70 efolds for successful inflation means that the multifield dynamics become 
negligible before perturbations on scales of observational relevance first cross the Hubble 
radius. Even if we push the observational window of interest back to iV* = 63 efolds before 
the end of inflation, rather than the usual assumption of iV* = 55 ± 5, we find that the 
model relaxes to effectively single-field dynamics prior to iV*. Hence the predictions from 
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Higgs inflation for observable quantities, such as the spectral index of the power spectrum 
of primordial perturbations, reduce to their usual single-field form. Moreover, the absence 
of multifield effects for times later than N* means that this model should produce negligible 
non-Gaussianities during inflation, in contrast to the broader family of models studied in 
[23]. 

The methods we introduce here may be applied to any multifield model with nonminimal 
couplings and an SO (AT) symmetry among the fields in field space. The conclusion therefore 
appears robust that such highly symmetric models should behave effectively as single-field 
models, at least within the observational window of interest between iV* = 63 and iV* = 40 
efolds before the end of inflation. Of course, multfield effects could become important in 
such models at the end of inflation, during epochs such as preheating. Such processes remain 
under study. 



APPENDIX A: ANGULAR EVOLUTION OF THE FIELD 



For completeness, let us integrate the angular equation of motion, Eq. (53), for all values 
of 7 (in the slow roll regime of the radial field). This yields 



7(t) (VXM pl +J2^ + XM^ 



7o 



cxp 



3A Afpit 



2s c 



(74) 



In the two limits, 70 <C y/\M p \/ '^/2£ and 70 3> V~\M p i/y/2£, we may solve Eq. (74) and 



recover the forms of ^(t) presented in Eqs. (55) and (57) 
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